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Abstract 

A three-component nonlinear Schrodinger-type model which describes spinor Bose- 
Einstein condensate ( BEC) is considered. This model is integrable by the inverse 
scattering method and using Zakharov-Shabat dressing method we obtain three types 
of soliton solutions. The multi-component nonlinear Schrodinger type models related 
to symmetric spaces C.I ~ Sp(4)/U(2) is studied. 

1 Introduction 

The dynamics of spinor BEC is described by a three-component Gross-Pitaevskii (GP) 
system of equations. In the one-dimensional approximation the GP system goes into the 
following nonlinear Schrodinger (MNLS) equation in (ID) x-space [1]: 

idt^i + dl^i + 2(|«>i|2 + 2|«>o|2)$i + 2$li$g = 0, 

idt^o + dl^o + 2(|^>_i|2 + |$oP + + 2$^^>i$_i = 0, (1) 

idt^-i + dl^-i + 2(|«>_i|2 + 2|$o|^)^-i + 2$I«>g = 0. 

We consider BEC of alcali atoms in the F = \ hyperfine state, elongated in x direction 
and confined in the transverse directions y, z by purely optical means. Thus the assembly 
of atoms in the F = \ hyperfine state can be described by a normalized spinor wave vector 

t) = ($i(x, t), ^>o(x, t), ^>_i(x, t)f (2) 

whose components are labelled by the values of rap = 1,0,-1. The above model is 
integrable by means of inverse scattering transform method [Ij. It also allows an exact 
description of the dynamics and interaction of bright solitons with spin degrees of freedom. 
Matter-wave solitons are expected to be useful in atom laser, atom interferometry and 
coherent atom transport. It could contribute to the realization of quantum information 
processing or computation, as a part of new field of atom optics. Lax pairs and geometric 
interpretation of the model ([I]) are given in [1]. Darboux transformation for this special 
integrable model is developed in [5]. The aim of present paper is to show that the system 
([1]) is related to symmetric space C.I ~ Sp(4)/U(2) (in the Cartan classification [8]) with 
canonical Z2-reduction and has natural Lie algebraic interpretation. The model allows 
also a special class of soliton solutions. We will show that they can be obtained by a 
suitable modification of the generalization of the so-called "dressing method" , proposed 
in [9]. 
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2 Solving the model: algebraic and analytic aspects. 

The model ([I]) belongs to the class of multicomponent NLS equations that can be solved 
by the inverse scattering method [6] . It is a particular case of the MNLS related to the 
C.I type symmetric space Sp(4)/U(2) [4j. These MNLS systems allow Lax representation 
with the generalized Zakharov-Shabat system as the Lax operator: 

Lip{x, t, ^) = ^-^ + (Qi^i - ^J)i^{x, t, A) = 0. (3) 

where J and Q{x, t) are 4x4 matrices: J = diag (1, 1, —1, —1) and Q{x, t) is a block-off- 
diagonal matrix: 

^ \pix,t) y' ' ' \<i>-i{x,t) -^o{x,t) 
p(xt)=( '^o^^'*) (4) 

Solving the direct and the inverse scattering problem for L uses the Jost solutions (j) = 
{4>~^, 4>~) and ip = {ip~ of ([3]) which are defined by, see [H] and the references therein: 



lim (/.(x,t,A)e^^'^^ = 1, lim ^^(2;, t, A)e*^-'^' = 1 (5) 



These definitions are compatible with the class of smooth potentials Q{x,t) vanishing 
sufficiently rapidly at a; — > ±00. It can be shown that (p'^ and ip'^ (resp. (p^ and ip~) 
composed by 4 rows and 2 columns are analytic in the upper (resp. lower) half plane of 
A. The scattering matrix associated to ([3]) is defined as 

T(t.A) = «,(.,t.A))-v(..t.A) = (;:(j;^^> -;:,<;-^-;' 

where a^{t, A) and b^{t, A) are 2x2 block matrices. The blocks a^, 6^, and satisfy 
a number of relations [111 [T2] ; for example 

a+(A)c-(A) + r(A)d+(A) = 1, a+(A)d'(A) - 6~(A)c+(A) = 0, (7) 

etc. The fundamental analytic solutions (FAS) x^{x:t,X) of L{X) are analytic functions 
of A for Im A ^ and are related to the Jost solutions by: 

X^ix, t, A) = ^{x, t, X)Sj{t, A) = i;{x, t, X)Tj{t, A). (8) 

Here S^, upper- and lower- block-triangular matrices: 



c+u (^d-{t,X)\ ^-(. ( c-it,X) 

^) - ^ c+(t, X))' ''J - V -d+{t, A) U 

^J^^^^)-[0 a-(t,A) )'^J^^'^^-[b+it,X)l 
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satisfying Tf{t,X)Sj{t,X) = T{t,X) and can be viewed as the factors of a generalized 
Gauss decompositions of T{t, A) [13]. If Q{x, t) evolves according to ([1]) then the scattering 
matrix and its elements satisfy the following linear evolution equations: 

i^ + 2AV(t,A) = 0, ^^ = 0, (9) 

so the block-matrices a^(A) can be considered as generating functionals of the integrals 
of motion. The fact that all 4 matrix elements of a"'' (A) for A S C+ (resp. of a~(A) for 
A G C_) generate integrals of motion reflect the superintegrability of the model and are 
due to the degeneracy of the dispersion law of ([I]). 

The system ([T]) can be written in a Hamiltonian form by introducing the Poisson 
brackets: 

{qj{x),Pk{y)} = 2i5kj5{x - y), {qi2{x),pi2{y)} = i^ix - y), (10) 
and the Hamiltonian H = -ffkin + -f^int^ 

i^kin = r dx + i + 

™ J~oo \ dx dx 2 \ dx dx dx dx J J ^ 

oo 



/oo 
-oo 

/oo 
dx (1^0^-1 + ^1^1 
-oo 



As mentioned above, one can use any of the matrix elements of a^(A) as generating 
functional of integrals of motion of our model. Generically such integrals would have 
non-local densities and will not be in involution. 

The classical i?-matrix approach [6l H] is an effective method to determine the generat- 
ing functionals of local integrals of motion which are in involution. From it there follows 
that such integrals are generated by expanding Inm^(A) over the inverse powers of A, see 
|13j . Here m^(A) are the principal minors of T(X); in our case 

m^{X) = a|i(A), rn^iX) = deta+(A), 

ni^ (X) = a22{X) , 7712 (X) = det a" (X) . i^"^) 
If we consider 

oo 

ln7n+{X)=^X-'^li'\ 
s=l 

then one can prove that the densities of /i*^'' are local in Q{x,t). The fact that |13j : 
{m±(A),m±(/x)} = 0, forA:,j = l,2, 

and for all A,^ G C± allow one to conclude that {I^\lp'^^ = for all k,j = 1,2 and 
s,p > 1. 

(2) (2) 

In particular, the Hamiltonian of our model is proportional to /g , i.e. H = 8ili^ . 
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3 Soliton solutions for the spinor BEC: The so(5) connection 

The soliton solutions of the sp(4) MNLS ([1]) were derived by using the dressing method [3]. 
They can be considered as particular cases of the soliton solutions of the generic MNLS 
eqs., derived through the matrix version of the Gel'fand-Levitan-Marchenko equation, see 
[21 [Ills]. Here we extend further these results and combining the ideas of [3l[l5] we specify 
three types of solitons for the model ([1]). 

We start our analysis with the well-known isomorphism between the algebras sp(4, C) 
and so(5, C) [8]. Since the Lax representation is of pure algebraical nature it is natural 
to expect that our model ([T]) can be treated also by an equivalent Lax operator L' whose 
potential Q'{x,t) and J' take values in so(5). A consequence of the above-mentioned 
isomorphism is that the typical representation of sp(4) used above is equivalent to the 
spinor representation of so(5). 

So we first remind some of our results in [10^ [TT| I12j . where we have constructed 
the fundamental analytic solutions, the dressing factors, the soliton solutions etc. for a 
class of Lax operators (including L'), related to the simple Lie algebra g, in the typical 
representation of g. So we first have to specify (if necessary) g ~ so(5) and then reformulate 
the corresponding results for the spinor representation of so(5). 

The main goal of the dressing method is, starting from a solution Xoix,t, A) of Lq{X) 
with potential Q(^Q'^{x,t) to construct a new singular solution xfix,t,X) with singularities 
located at prescribed positions Xf; the reduction p = used in eq. (jH) ensures that 
Aj~ = (A]*")*. The new solutions xf{x, t, A) will correspond to a potential Q(i){x, t) of L{X) 
([3]) with two discrete eigenvalues Xf. It is related to the regular one by a dressing factor 
u{x, A) 



xf{x,t,X)=u{x,X)xQ{x,t,X)u_^{X). u_(A) = lim u{x,X) (13) 



Note that n_(A) is a diagonal matrix. The dressing factor u{x, A) must satisfy the equation 

.du 

dx 



du 

i— + {x)u - uQ(^o) (x) - X[J, u{x, A)] = 0, (14) 



and the normalization condition limA-+oo u{x, A) = U. Besides xf{x, A), i = 0, 1 and u{x, A) 
must belong to the corresponding Lie group 5p(4, C); in addition u{x,X) by construction 
has poles and/or zeroes at A^. 

The construction of u{x, A) is based on an appropriate anzats specifying explicitly the 
form of its A-dependence. 

n(x,A)= ll + (ci(A)-l)Pi(x,t) + f^-l)Pi(x,t), 



ci(A) 

ci(A) = ^f^, (15) 

where the projectors Pi{x,t) and Pi{x,t) are are of rank 1 and are related by Pi{x) = 
SP[{x)S-^. They must satisfy Pi{x,t)Pi{x,t) = Pi{x,t)Pi{x,t) =0. By S we have 
denoted the special matrix which enters in the definition of the orthogonal algebra, i.e. 
X G so(5) if X + SX^ = 0. In the typical representation of so{5) we have S = 
Tl=i{-^f^^Ek,^-k where (^y) 

km — ^ik^jm- The Construction of Pi{x,t) and Pi{x,t) 
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using the 'polarization' vectors is outlined in [llj and we skip it. The new potential is 
obtained from 

g(i) {x, t) - Q(o) {x, t) = (A+ - A^) [ J, Pi {x, t) - Pi (x, t)] . 

Here we show that the A-dependence of u{x, A) may depend [10] on the choice of the 
representation of so(5,C) ~ sp(4, C)). For so(5) it was shown |IO l \TT l [T2 | [T5] that there 
are three types of solitons: 



• the first type of soliton solutions are generated by dressing factors of the form (jlSp . 
For generic choice of the polarization vectors Pi{x,t) — Pi{x,t) £ so(5). 

• the second type of soliton solutions are generated analogously with dressing factor 
(fT5]) . but due to a specific choice of the polarization vectors Pi{x,t) — Pi{x,t) E 
so(3) C so(5). 

• the third type of soliton solutions are generated again by (jlSp but now the corre- 
sponding projectors Pi{x,t) and Pi{x,t) have rank 2. 

Each of these types of soliton solutions have their counterpart relevant to our model on 
sp(4). To the first type of soliton solutions there correspond dressing factor and potential 
Q(^i-^{x,t) of the form[5]: 

u{x,X) = ^y ci{X)tti{x) + J— - 7fi(a:), 

VCi(A) 

Q{l){x,t) - Q{0){x,t) = -[J,TTi{x,t) - 7fi(x,t)], (16) 

where vri(x) and vfi(x) are rank 2 projectors, such that 

7fi(x)7ri(x) = 7ri(x), 7fi(x)=0, 7fi(x) + 7ri(x) = U. (17) 

This last property ensures the non-degeneracy of u{x,X). Note that now the dressing 
factor is not a rational function of A but for the dressed FAS x{x, A) eq. (jl3p we get: 

xf{x,t,X) = ^vri(x,t) + — ^7fi(x,t)^ Xo{x,t,X) (vrf + ci(A)7ff) , 

7r7 = lim 7ri(x,t), (18) 

i.e., the fractional powers of ci(A) disappear. 

The second type of solitons with rank 2 projector Pi (x) after recalculating to the spinor 
representation formally keeps the same form (llSp with Pi(x) replaced by ^i(x) which has 
rank 1 but generically is not a projector, see [Sj. 

The third type of solitons is similar to the second one but with additional constraints 
on the factor ^i(x) so that Ai{x) — Ai{x) G sp{2) C sp(4). 

Consider the purely solitonic case when Q(o) = 0- From now on we introduce the 
following notations A^ = /ii it ii^i and 

A = -2i((A+)2 - (Ar)2)t - i(A+ - Ar)x, 

B = -2{{Xt f + (Ar)2)t - (A+ + X^)x. (19) 
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Here A{x,t) and B{x,t) are x and t dependent real valued functions. Making use of 
the explicit form of the projectors P±i{x) valid for the typical representations of B2 we 
obtain [11]: 

$(1) (x, t) = 7 \^ ^ {no,imo.2e^ + Uq -^rriQ le'^) e'^ (20) 
<l>(0)(x,t) = ^^2/M—^ (no,imo,3e^ - n,-,m,-,e-^) e^^ (21) 

$(_i)(x,t) = -li^i-ZAl (no,imo,2e^ + no,2moje-^) e'^ (22) 

where the denominator in the above formula is given by: 

{m\n) = mo,ino,i(e^^) + mo_inoj(e"^^) 

+ "^o,2no,2 + "Zo,2"-o,2 + "T-o.s^-cs- (23) 

and mo,fc, no,fc are the components of the polarization vectors. 

Choosing appropriately the elements of the polarization vectors |no) and |mo), one can 
show that the conjecture that the Zakharov-Shabat dressing procedure and the Gel'fand- 
Levitan Marchenko formalism lead to comparable soliton solutions is true. It is not a 
problem to multiply the polarization vectors |no) and \mo) by an appropriate scalar and 
thus to adjust the two solutions. Such a multiplication easily goes through the whole 
scheme outlined above. The involution Qj-j^^ = that the potential of the Lax operator 
([3]) associated with the system ([T]) is subject to results in the following relations between 
the elements of the "polarization" vectors |no) and (mo|, namely UQ^k = "^ofc- Utilizing 
the above and a proper change of field components, we can relate the solution 

Cte^ + a2CV2 det{Ct}e- ^ 
+ Ty + |det{C7}|2e-2^ 

where W = (2|ci2p + jcip + |c2p) and the "polarization" matrix can be cast into the form 



C=r^2 (25) 

V C2 -Cl2 J 

In the special case when W = 1 and det{C} = we obtain 
cosh A 

Thus we confirm the result obtained in[lj, aquaired with the help of GLM formalism 
and the solution (j20p . derived within the generalized Zakharov-Shabat dressing procedure, 
provided we make sure that the extra condition on the vector \m): 

- 2mo,imo_i + 2mo,2"io,2 = ("lo.s)^ , (27) 
and analogous one for |n) holds true. Setting 

mo,i = 1, iriQ i = -(c*2)^ - c*C2, 

"10,2 = ici, 'm,o,2='^C2, mo^s = 'mo^s = -V2cl2 
we establish the equivalence between the two solutions. 
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4 Conclusions 

We have derived the sohton solutions of the three-component system of NLS type on the 
symmetric space Sp(4)/U(2) which is related to spinor Bose-Einstein condensate model 
(with F = 1). Furthermore, we have described briefly the Hamiltonian properties of the 
model and the integrals of motion. Using the classical r-matrix approach, we showed that 
the integrals of motion, that belong to the principal series are in involition. 

The reduction of the multi-component nonlinear Schrodinger (NLS) equations on sym- 
metric space C.I ~ Sp(2p)/U(p) for p = 2 is related to spinor model of Bose-Einstein 
condensate. Other interesting reductions of MNLS type equations were reported in [11] 
and a systematic study of the problem is on the way. Recently the authors of [16] develop 
a perturbation theory for bright solitons of the F = 1 integrable spinor BEC model. Both 
rank-one and rank-two soliton solutions are obtained using Riemann-Hilbert method and 
are compared with known results. 
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